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Abstract. We construct examples of 2-step Carnot groups related to quaternions and study 
their fine structure and geometric properties. This involves the Hamiltonian formalism, which 
is used to obtain explicit equations for geodesies and the computation of the number of 
geodesies joining two different points on these groups. We able to find the explicit lengths of 
geodesies. We present the fundamental solutions of the Heat and sub-Laplace equations for 
these anisotropic groups and obtain some estimates for them, which may be useful. 



1. Introduction 

This paper presents examples of Carnot groups and studies their fine structure, geometric 
properties and basic differential operators attached to them. A Carnot group is a connected 
and simply connected m-step nilpotent Lie group G whose Lie algebra Q decomposes into the 
direct sum of vector subspaces V\ ® V2 ffi . . . © V m satisfying the following relations: 

[Vi,V k ] = V k+1 , 1 < k < m, \Vi, V m ] = {0}. 

The simplest examples of this group are Euclidean space M. n , Heisenberg group HP and 
H(eisenberg)-type groups introduced by Kaplan |18| . The Carnot groups form a natural habi- 
tat for extensions of many of the objects studied in Euclidean space and find applications 
in the study of strongly pseudoconvex domains in complex analysis, semiclassical analysis of 
quantum mechanics, control theory, probability theory of degenerate diffusion process and oth- 
ers. The geometry of Carnot groups and differential operators related with them were studied 
extensively by many mathematicians, for instance, in pfl IH 151 1§1 ITU [TBI |2*U1 122] . 

We construct examples of 2-step Carnot groups, related to the multidimensional space of 
quaternion numbers. We will call these groups anisotropic quaternion groups and denote 
them by Q n . In ^2] the quaternion H-type groups were studied. The results of ^21 can be 
easily extended to the multidimensional quaternion space. The examples of the present paper 
contain the multidimensional quaternion H-type groups as a particular case. We construct 
the Hamiltonian function associated with the sub-Laplacian generated by left invariant vector 
fields. Solving the Hamiltonian system of differential equations we give exact solutions that 
describe the geodesies on the group. We study geodesic connectivity between any two points of 
the group. It is known that every point of a Riemannian manifold is connected to every other 
point in a sufficiently small neighborhood by one single, unique geodesic. But in this case, 
there will be points arbitrarily near a point which are connected to this point by an infinite 
number of geodesies. Since we are working on a group, we may simply assume that the point 
is the origin O = (0,0). We prove the following results: 

2000 Mathematics Subject Classification. 53C17, 53C22, 35H20. 

Key words and phrases. Quaternions, Carnot-Caratheodory metric, nilpotent Lie groups, Hamiltonian for- 
malism, Green function. 

This work was supported by Projects FONDECYT (Chile) # 7050181, #1040333, and by grant of the 
University of Bergen. 



2 



DER-CHEN CHANG, IRINA MARKINA 



(1) . If P = (x, 0) = (xi, . . . ,x n , 0) with x / 0, then there is only one geodesic connecting 
the origin O and the point P; 

(2) . If P = (xi, . . . , x n , z) with X[ 0, I = 1, . . . , n and 2^0, t/ien i/iere are finitely many 
geodesies connecting the origin O and the point P; 

(3) . If P = (xi, . . . , x n , z) 7^ /or / = 1, . . . ,p — 1 and x; = for I = p, . . . ,n, z / 0, 
then there are countably infinitely many geodesies connecting the point O and the point P; 

(4) . If P = (0, z) with z ^ 0, then there are uncountably infinitely many geodesies connecting 
the point O and the point P. 

We will discuss basic properties of geodesies in sections 2 to 4. Then we will prove connec- 
tivity theorems in section 5 (see Theorems 15.11 I5~2l 15.41 and 15 .6|) . Furthermore, parametric 
equations and arc lengths of all these geodesies are calculated explicitly. 

We also consider complex geodesies and find a relation between the complex action function 
and the Carnot-Caratheodory metric. The complex action function allows us to deduce the 
transport equation and its solution: volume element. The fundamental solution of the Heat 
equation is given in terms of the complex action function and volume element. More precisely, 
the heat kernel at the origin is given by 



C f -f 



where 
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f(y,w,r) = -i^2T m w m + ^^p|r|;coth(|r| 

m 1=1 

is the modified complex action and 



V(r) = n 



l=1 sinh 2 (|r|/) 

is the volume element. (See Theorem I7.1JI . Integrating the fundamental solution of the Heat 
equation with respect to the time variable, we obtain the Green function for the sub-Laplacian 
(see Theorem I7.2JI : 

_ 2^ (2^"+3 r V[T + l£ ~ z) 
{ ' j " (2n + l)! J R3 pn+2 {T + i£2 ) dT - 
The last section is devoted to some estimates of fundamental solutions, which may be useful. 
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2. Definitions 

A quaternion is a mathematical concept (re) introduced by William Rowan Hamilton from 
Ireland in 1843 j^j- (It has been said that when Hamilton discovered the quaternions, they 
stayed discovered). The idea captured the popular imagination for a time because it involved 
relatively simple calculations that abandon the commutative law, one of the basic rules of 
arithmetic. Specifically, a quaternion is a non-commutative extension of the complex numbers. 
As a vector space over the real numbers, the quaternions have dimension 4, whereas the 
complex numbers have dimension 2. While the complex numbers are obtained by adding the 
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element i to the real numbers which satisfies i 2 = —1, the quaternions are obtained by adding 
the elements i, j, and k to the real numbers which satisfy the following relations 

i 2 = j 2 = k 2 = ijk = -1. 

Unlike real or complex numbers, multiplication of quaternions is not commutative, e. g., 

(2.1) ij = -ji = k, jk = -kj = i, ki = -ik = j. 

The quaternions are an example of a division ring, an algebraic structure similar to a field 
except for commutativity of multiplication. In particular, multiplication is still associative and 
every non-zero element has a unique inverse. 

The quaternions can be written as a combination of a scalar and a vector in analogy with 
the complex numbers being representable as a sum of real and imaginary parts, a ■ 1 + b ■ i. For 
a quaternion h = a + bi + cj + dk we call a scalar a the real part and the 3-dimensional vector 
u = bi + cj + dk is called the imaginary part of h and is a pure quaternion. In M 4 , the basis of 
quaternion numbers can be given by real matrices 



U = 



10 
10 
10 
1 



Mi 



We have 
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aU + bMi + cM 2 + dM 3 . 



Similarly to complex numbers, vectors, and matrices, the addition of two quaternions is equiv- 
alent to summing up the coefficients. Set h = a + u, and q = t + xi + yj + zk = t + v. 
Then 

h + q = (a + t) + (u + v) = (a + t) + (b + x)i + (c + y)j + (d + z)k. 

Addition satisfies all the commutation and association rules of real and complex numbers. The 
quaternion multiplication (the Grassmanian product) is defined by 

hq = (at — u • v) + (av + tu + u x v), 

where u • v is the scalar product and u x v is the vector product of u and v, both in M 3 . The 
multiplication is not commutative because of the non-commutative vector product. The non- 
commutativity of multiplication has some unexpected consequences, e.g., polynomial equations 
over the quaternions may have more distinct solutions then the degree of a polynomial. The 
equation h? + 1 = 0, for instance, has infinitely many quaternion solutions h = a + bi + cj + dk 
with b 2 + c 2 + d 2 = 1. The conjugate of a quaternion h = a + bi + cj + dk, is defined as 
a — bi — cj — dk and the absolute value of h is defined as \h\ 



h* = a — bi — cj — dk and the absolute value of h is defined as \h\ = Vhh* = \Ja? + b 2 + c 2 + d 2 . 

Let us denote the space of quaternions by 7i. We consider n-tuples of quaternions: H = 
(hi, I12, ■ ■ ■ , h n ), hi e H, i = 1, . . . , n. We may define addition between two of them in an 
obvious way: 

H + Q = (hi + qi, h 2 + q 2 , ■ ■ ■ , K + q n ) 
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for H = {h\, . . . ,h n ) and Q = (qi, . . . , q n ). Multiplication by scalar is defined as aH = 
(ah\, ah 2 , ■ ■ ■ , ah n ), where a may be real or complex number. Therefore, the n-dimensional 
quaternions, Ji n is a vector space. The norm is defined by 



n 

\h\w = (En 2 ) 



1/2 



In this article we will construct 2-step Carnot groups related to the multidimensional quater- 
nion numbers. We will call these groups anisotropic quaternion groups Q n . To explain precisely 
the idea of their construction we first introduce the notion of H-type Carnot groups and then, 
making some modifications, we arrive at our main example. 

H-type homogeneous groups are simply connected 2-step Lie groups G whose algebras Q are 
graded and carry an inner product such that 

(i) Q is the orthogonal direct sum of the generating subspace V\ and the center V 2 : 

g = v l ® V 2 , V 2 = [V u Vx], [V u v 2 ] = 0, 

(ii) the homomorphisms Jz '■ V\ — ► Vi, Z G V 2 defined by 

(JzX,X') = (Z,[X,X']), X,X'€V! 

satisfy the equation 

Jf = -\Z\ 2 I, Z G V 2 . 
Here (•,•) is a positive definite non-degenerating quadratic form on Q, [-,-] is a commutator 
and I is the identity. The group is generated from its algebra by exponentiation. 

To construct the multidimensional quaternion H-type group we take the space of quaternions 
7i n as V\ and generate the center V 2 . We consider the n-dimensional imaginary quaternions 

Z x = (aqi, . . . ,aii), Z 2 = (a 2 j, . . . , a 2 j), Z 3 = (a 3 k, . . . , a 3 k) 
with positive constants a m . They have the following representation as real matrices 4n x 4n: 

a m M m 



M m = 







a m Mr 



where there are n blocks on the diagonal of each matrix M m , m = 1,2,3. The matrices M m , 
m = 1,2,3, are the matrices associated to the homomorphisms Jz- 

Now we extend the construction, introducing anisotropy to this very symmetric setting. We 
take an arbitrary n-dimensional imaginary quaternions 



Z x = (ani, . 

with a m i > for all m = 
are following: 



Mi = 



auMi 



. . ,ai„i), Z 2 
1, 2, 3 and I = 





= («2lj, • 

1, . . . ,n. 



■ ,a 2n 'i), Z3 = (a 3 ik, . . . ,a 3n k), 

The representation as real matrices 4n x 4n 



M 2 = 



a ln Mi 



a 2 \M-i 







a 2n M 2 



M« = 







where there are n blocks on the diagonal of each matrix M m , m = 1,2,3. We construct 
and make the principal calculations for the anisotropic quaternion group Q n with center V 2 of 
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topological dimension 3. For the other examples see Remark 12.11 The corresponding algebra 
Q n is the two-step algebra Vi©V2. The topological dimensions of the group is dimQ™ = 4n+3. 
The homogeneous dimension defined by the formula v = dim V\ + 2 dim V 2 plays an important 
role in analysis on homogeneous groups. We see that the homogeneous dimension is always 
greater than the topological dimension and in our case equals v(Q n ) = 4n + 6. An isotropic 
case based on the one dimensional space of quaternions was studied in ^2] an d the anisotropic 
Heisenberg group was considered in 

We set the standard orthonormal systems {X^i} G K 4n , k = 1,2,3,4, I = l,...,n and 
{Z\, Z21 Z 3 } G M 3 . We reserve the following indexes: I = 1, . . . n denotes the coordinate index 
in (h\, . . . , h n ); k = 1, 2, 3, 4 denotes the index of coordinates inside of each quaternion hi and 
m = 1,2,3 is related to the coordinate index in V2 or the index of the matrices M m . The 
matrices M m transform the basis vectors in the following way 
(2.2) 

MiXu = -a u X 2 i, MiX 2l = auXn, MiX 3l = -a n X ih = a u X 3l , 

M-2X11 = (121X41, WI2X21 = 0,21X31, M2X3/ = —CI21X21, M.2X41 = —a,2iXu, 
M 3 X U = a 3l X 3l , M 3 X 2 i = -a 3l X 4l , M 3 X 3l = -a 3l Xu, M 3 X 4/ = a 3l X 2l . 

We use the normal coordinates 

q = (x,Z) = (xii,X21,X 3 i,X 4 i, . . . , Xi n ,X2n-,X 3n , X 4n , Z X , Z 2 , Z 3 ) 

for the elements 

exp ^ XkiX ki + ^2 z mZ m ^ G Q n - 
k,l m 
The Baker-Campbell-Hausdorff formula 

exp(X + Z) exp(X' + Z') = exp(X + X' , Z + Z' + ^[X, X']), 

for X, X' G V\, Z, Z' G V2 defines the multiplication law on Q n . Precisely, we have 
L q (q') = L^ z ) (x, z) = (x, z) o (x 1 , z ) 

= (x + x',zi + z[ + ~(Mix,x / ),Z2 +z' 2 + ^(M 2 x,x'),z 3 + z 3 + ^(M 3 x,x')), 

for q = (x,z) and q' = (x',z'), where (M m x,x') is the usual scalar product of the vector 
M m x G M 4n by x' G M 4n . The multiplication "o" defines the left translation L q of q' = (x' , z') 
by the element q = (x, z) on the group Q n . 

We associate the Lie algebra Q n of the group Q n with the set of all left invariant vector fields 
of the tangent bundle TQ n . The tangent bundle contains a natural subbundle TQ n consisting 
of "horizontal" vectors. We call TQ n the horizontal bundle. The horizontal bundle is spanned 
by the left-invariant vector fields X\\(x, z), . . . , X 4n (x, z) with Xki(0, 0) = Xki, k = 1, . . . , 4, 
/ = 1, ... ,71, (see for example, jZ], [Hj, |2H])- In coordinates of the standard Euclidean basis 
gf-, these vector fields are expressed as 

8 1/ 

Xn(x, z) =- h - +aux 2 i 

ox\i 2 V 



X 2 i(x,z) 



d 1 



(2.3) 



dx 2 i 2 



d 1 

X 3i (x,z) =- h -( +a lt x A i 

dx 3 i 2 

d 1 

X<u(x,z) =- h - ( - aux 3l 
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- a 2 iX4i-— 
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- a 3 ix 3 i — 
dz 3 
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- a 2 ix 3 iT— 
0Z2 


+ a 3 ix 4 i — 
oz 3 
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dz\ 


+ a 2 lX2lT— 
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+ a 3 ixn — 
oz 3 
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- a 3 ix 2 i^— 
dz 3 
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for Z = 1, ... j 7i. The left invariant vector fields Z m (x, z) with Z m (0, 0) = Z m , m = 1, 2, 3, are 
simply the vector fields 

(2.4) z m {x,z) = ^-. 

oz m 

We write simply X^i and Z m instead of Xki(x,z) and Z m (x,z), if no confusion may arise. 
Note that if we fix m to be equal only to 1,2 or 3 then the vector fields (|2.3fl are reduced 
to the anisotropic vector fields of anisotropic H 2n Heisenberg group and the group Q n is 
isomorphic to anisotropic H 2n Heisenberg group, considered in |Sj. We also use the notation 
Xi = (X\i, . . . , X41), I = 1, . . . , n. We call the next vector 

X = (X 11: ...,X 4n ) = (V x + - ^M m x— - ), 

z oz m 

m=l 

where V x = (g§^, ■ ■ ■ , q® a ) the horizontal gradient. Any vector field Y belonging to TQ n is 
called the horizontal vector field. In particular, the horizontal gradient X is a horizontal vector 
field, that justifies the name "horizontal" gradient. 
The commutation relations are as follows 

[Xu,X 2 i] = -auZi, [Xu,X 3 i] = a 3 iZ 3 , [Xu,X 4 i] = a 2 iZ 2 , 

[X2i,X 3 t] = a 2 iZ 2 , [X 2 i,X4i] = —a 3 iZ s , [X 3 i,X 4 i] = —auZi, 
for any I = 1, . . . , n. 

A basis of one- forms dual to X^, Z m , is given by dxki,$ m , with 

(2.5) -dm = dz m - - (M m x, dx) . 

Since the interior product "& m (X).i) vanishes for all m = 1,2,3, k = 1, . . . , 4, I = 1, . . . , n we 
have the product i? m (Y) vanishing on all horizontal vector fields Y . 

Remark 2.1. If we formally put an = 0, / = 1, . . . ,n, then we obtain another example of a 
quaternion anisotropic group with 2-dimensional center. The case an = a 2 i = 0, I = 1, . . . ,n, 
corresponds to an anisotropic group with 1-dimensional center. 

3. Horizontal curves and their geometric characteristics 

Summarizing the results of the previous section we can say that Q n is a space of (4n + 3) - 
tuples of real numbers R 4n+3 where the commutative group operation "+" is replaced by 
the noncommutative law "o". Respectively, the left translation L x {x') = x + x' (that in the 
commutative case coincides with the right translation) is substituted by the left translation 
L q (q') = q o q' . The corresponding Lie algebras are fundamentally different. The constant 
vector fields Xi = t£-, i = l,...,4n + 3, of Euclidean space are replaced by the vector 
fields (|2.3|) and (|2.4j) . Moreover, since the group of vector fields (|2.4|) is completely generated 
by the group (|2.3j) by means of commutation relations, the geometry of the group Q n is defined 
by the horizontal bundle TQ n . The velocity and the distance should respect the horizontal 
bundle TQ n . Since [X^i , Xk a i b ] ^ TQ n the horizontal bundle is not integrable, i.e., there is no 
surface locally tangent to it [Q. As we say, the geometry is defined by the horizontal bundle, 
so it is sufficient to define the Riemanian metric only on the horizontal bundle TQ n of tangent 
bundle of Q n . These kinds of manifolds have aquired the name subRiemanian manifolds. The 
definitions and basic notations of subRiemanian geometry can be found in, e.g., |24| . 

A continuous map c(s) : [0, 1] — ► Q n is called a curve. We say that a curve c(s) is horizontal 
if its tangent vector c(s) (if it exists) belongs to TQ n at each point c(s). In other words, there 
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are (measurable) functions a k i(s) such that c(s) = Ylk i a ki( s )X k i(c(s)). We present some 
simple propositions that describe the geometry of the group Q n . 

Proposition 3.1. A curve c(s) = (x(s), z(s)) is horizontal if and only if 

1 

(3.1) z m = -(M m x,x), m = 1,2,3, 

Where X = (in, ±21, X3I, ±41, • • • , ±ln, £2n, ^3n, ±4n)- 

Proof. We can write the tangent vector c(s) in the form 

d d 
c(s) =(x(s), z(s)) = x k i(s)- h i m (s x 



OX hi oz m 
k,l m m 



1 d \ / 1 \ 5 

m m m 



It is clear that c(s) is horizontal if and only if the coefficients in front of gj— , m = 1,2,3, 
vanish. This proves Proposition 13.11 □ 

Corollary 3.2. If a curve c(s) = (x(s), z(s)) is horizontal, then 

c(s) = (x(s),x\ = ^2x kl (s)X kl . 

k,l 

It is easy to see the following statement. 

Proposition 3.3. The left translation L q of a horizontal curve c(s) = (x(s), z(s)) is a hori- 
zontal curve c(s) = L q (c(s)) with the velocity 

(3.2) 5( s ) = (L,),c(s) =J2^kl{s)X kl (c(s)) = (x(s), X(c(s))) . 

k,l 

Proposition 3.4. The acceleration vector c(s) of a horizontal curve c(s) is horizontal. 

Proof. Let c(s) be a horizontal curve. Then c(s) G TQ r ^ s y Let us show that c(s) E ^~Qc( s )- 
Differentiating equalities (|3.1|l of the horizontality condition and making use of (Wl m x, x) = 
for m = 1, 2, 3 and any x £ M 4n , we deduce that 



i((M m a;(s),s(s)) + (M m x{s),x{s))^j = -(M m x{s),x{s)) 
for m = 1,2, 3. Then the acceleration vector along c(s) is 

c(s) =(£(*), V*) + (ij(s),V z ) = (i(s),V x . + -^M m x(s)-— ) 

m m 

+ E (^( s ) " ^(^),M m x( S )))^ = (x(s),X(c(s))). 

m m 

This means that the vector c(s) is horizontal. The proposition is proved. □ 



8 der-chen chang, irina markina 

4. Hamiltonian formalism 

In this section we study the geometry of the anisotropic quaternion group Q n making use 
of the Hamiltonian formalism. The geometry of the group is induced by the sub-Laplacian 
Ao = Ylkl-^-kl- Operators of such type are studied, for instance, in [1JEI]. Since the vector 
fields X^i satisfying the Chow's condition, by a theorem of Hormander ^7j, the operator Aq 
is hypoelliptic. Explicitly, the sub-Laplacian has the form: 

3 n r,o 3 



where V x 



n 4 

Y.Y. 

i=i fc=i 



X 



hi 



A 



m=l 



^2 a2 ml\ x l\ 2 ) 
1=1 



q2 ° Q 



ti.i'i 



and \xi\ 2 = Y, 



m=l 
2 



dz r , 



To present the 



Hamiltonian function we introduce the formal variables £ = (£n, . . . , 6in) with £ 



kl 



h, O3) with 9 m = tM-, m = 1, 2, 3. The associated with sub-Laplacian Ao Hamiltonian 



function H(£, 9, x, z) is following 



(4.1) 



H(£,9,x,z) = \H\ 2 + 



-.3 n 



2 I |2 
, a ml\ x l\ 



+ 



m=l 1=1 



where |£| 2 = Yli k ^tv anc ^ diagonal blocks of the matrix ]P 



6 m M m )x,£), 

are of the form 



m=l 








— 03l93 


—02/6*2 


an9i 





— 02l92 


03163 


0.3163 







a>il9i 


02162 


— a 3l63 


—o\i9\ 






We use the following notation: 



ml 



©2 

convenience: 



Zt=i6 2 m A 2 m , a ndM = Z 
2 



Em=l 6 



2 a 2 



ra=l rn u 'ml > 





x 6* m M r 



a mrM 



We also introduce some different metrics for 



In this notation we get Y^i=i o 2 
this notation 

(4.2) H(£,e,x,z) 



ml\ X l\ 



{B x,x) = (Bx,Bx), where B is a diagonal matrix. 
. The Hamiltonian function takes a new form in 



1^ 



lei 2 + \ E Ml 



4 + (Mx,e) = ici 2 + ^(© 2 x,x) + (Mx,e), 



and the corresponding Hamiltonian system obtains the form 

= ^ = 2£ + Mx 



(4.3) 



---7/ 



m 



i 



as 

<9Zm 



1,2,3. 



±0 2 x + M£ 



0. 



The solutions 7(5) = (x(s), z(s), £(s), 9(s)) of the system (|4.3j) are called bicharacteristics. 

Definition 4.1. Let Pi(xo,#o); P2( 2; > 2 £ Q n - A geodesic from Pi to P2 is the projection of 
a bicharacteristic 7(5), s G [0, r], onto the (x, z)-space, that satisfies the boundary conditions 

(x(0),z(0)) = (x ,z ), (x(t),z(t)) = (x,z). 
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The next properties of the matrices M m , m = 1, 2, 3, are obvious: 

(4.4) M 2 m = —U, m = 1,2, 3, where W is the unit (4x4)- matrix. 

MxM.2 = -M 2 M\ = M3, M2M3 = -M3M2 = Mi, 

(4.5) M3M1 = -MiM-z = M 2 . 

(4.6) -^m 1 = —Mm, where .M" 1 is the inverse matrix of M m , m = 1, 2, 3. 

(4.7) = —A4 m , where is the transposed matrix for A4 m , m = 1,2,3. 

(4.8) (M m x, x) = 0, m = 1, 2, 3, for any x e M 4 . 
As a corollary we obtain some useful formulas. 

Proposition 4.2. In the above-mentioned notations we have 

n 

(4.9) (M m x,Mx) = 9 m ^a 2 ml \xi\ 2 = d m \x\\ m for any j = 1,2,3. 



(4.10) M 2 = - d m A m = -® 2 > M 3 = -® 2 M, M 4 = 4 , M 5 = 4 M . . . 
Proof. 



m=l 



(M m x, Mx) = (M m x, ^2 6j M i x ) = E 0j(M m x, Mjx) 

3=1 3=1 

3 

= Qj ((ijlMjXi, djnMjXn) , (a m iM m xi,... ,a mn M m x n ) 
3=1 ' 

3 n 3 n 

= y^9jy^ajia m i(MjXi,M m xi) = / J 0j/ J ajia m l(-M m MjXi, xj) 

j=l 1=1 j=l 1=1 



n 



1=1 

by the properties (|Q|) . (|4~5|) . (|4~7|) . and (|4~8|) of matrix M m . 

To prove dHHj) we note that M m Mj = -MjM m by (|43|) and = -A^ by the prop- 
erty (|4.4|) for any j, m = 1, 2, 3. Then 

3 3 3 3 

M ' = E E ^mMjMm = £ ^M 2 ^ = - J] ftm^m = "© 2 . 
j'=l m=l m=l m=l 

The rest is obvious. □ 
Lemma 4.3. Any geodesic is a horizontal curve. 
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Proof. Let c(s) = (x(s), z(s)) be a geodesic. The system 1)4.3)1 implies 

(4-11) zm = ^rML + \ {M m x, 2£) = ^\x\ 2 Am + ^ (M m x, x) + ^ (M m x, 2£ - x) . 
Making use of the first line of the system 1)4.3)) , we write the last term of 1)4.11)) as 

(4.12) ~(M m z,2£-x) = -i(M m x,Mx) = -^\x\ 2 Am . 
Here we used the formula 1)4.9)1 . Combining (|4.11|) and (|4.12f) we deduce 

(4.13) i m = ^\x\ Am + (M m x,Z) = i(M m x,x), m = l,2,3. 

Therefore, c(s) is a horizontal curve by Proposition 13.11 □ 

Lemma 14.31 shows that the second equation of the system (|4.3|) is nothing more then the 
horizontality condition ()3.1)) . 

Let us try to solve the Hamiltonian system explicitly. The last equation in (|4.3|) shows that 
the function H(£,0,x, z) does not depend on z. We obtain that 6 m are constants which can 
be used as Lagrangian multipliers. Multiplying the first line of system (|4.3|) by M, we obtain 

(4.14) Mi = 2M£ - e 2 x. 

Expressing M£ from Q4.14JI and substituting it in the equation for £ from (|4.3|) . we get 
(4.15) 

We differentiate the first equation of 1)4.3)1 and substitute the £ from 1)4.15)1 . Finally, we deduce 

(4.16) x = 2£ + Mi = 2Mi. 

Let us solve the equation 1)4. 16|) . We substitute y(s) = x(s). The equation y(s) = 2My(s) 
has a solution y(s) = exp(2sM)y(0). Therefore, 

(4.17) x(s) = exp(2sM)x(0). 

Let us discuss the properties of the matrix exp(2sM). For simplicity of notation we write [B]i 
for Z-block of a block diagonal matrix B. 

Lemma 4.4. The exponent exp (2sM) is an antisymmetric block matrix that commutes with 
M and which blocks can be written in the form: 

(4.18) [exp(2 S M)] ; = cos(2 S |0|,)« + Sm ^ l) \M]i- 
Proof. We observe that 

n=0 k=0 y ' 1 14 fc=0 V ; 

^ {2s\9\ l f k+2 [Mjj ^ (2s\6\i) 4k+3 
^ (4fc + 2)! ^ (4fc + 3)! 

by ()4.1U)I . We conclude that the matrices M and exp (2sM) commute. Note that 

^ (2»i»ii)" fwr ...... 
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and 



^ (4fc + l)! ^ 



{2 S \e\i) 4k+3 

(4k + 3)\ 



sm(2s\6\i). 



k=0 v ' k=0 

With this notation, exp (2sM) is a block diagonal matrix with the blocks [exp (2sM)]^ 



cos(2s|0|/) S m(2s\e\i) -«^ci,wo t .k)i,, _«^k 



sin(2 S |^) -^sin(2 S |0|O 



•fsn(28\e\i) cos(2s|0|,) -2fgjf sin(2s|0|,) * sin(2 S |0|,) 

2j^sm(2 S |0|,) s in(2 S |0|,) cos(2 S |0|,) 2^ sin(2 S |0|,) 

^sin(2 S |0|,) -^fsin(2 S |^) sin(2 S |^) cos(2#|,) 



□ 



The group structure allows to restrict our considerations to the curves issuing from the 
origin. Hence, x(0) = 0. The equation (|4.17|) has the form 



(4.19) 



£,(*) = cos(2 s |0|,)W£,(O) + g^M [M ],£,(Q), l = l,...,n, 

Ml 



by (|4.18j> . Integrating from to s we get 

, /^ 1 - cos(2s|6»|i) r , ^ . . sin(2s|6»| z ), < , N 
(4.20) Xl (s) = 2|g y 1 U; [M],£,(0) + _^J-i^£,(0), Z = 1, . . . , n. 

Let us describe the z-components of a geodesic curve. If a curve is geodesic, then it is horizontal 
by Lemma ESI and we have 

■ n ^tv/t r ^ 'M^ yy cos(2a|gjj)(l - cos(2s|flLJ) , . . , 

^m(s) = -(M m x(s),x(s)) =Z^{ Tmp ([M m ],[M],x,(0),a:,(0)) 



+ 



4|*l? 

sin(2a|g|i)(l-cos(2s|g|Q) 
4|0|? 



([M m ],[M],£,(0),[M],£,(0)) 



+ sm(2^|Ocos ( 2 S |^ ) ([Mmhx;(o);xKo)) 



+ 



4|<9|i 
sin 2 (2#l,) 
4|0| 2 1 



[M m ],£,(0),[M],£,(0)) , 



for m = 1,2, 3 by (|4~H)|) and The properties l|4TH|l and l)OJl imply 

([M m ],£,(0),£,(0)) = ([M OT ],[M],£,(0),[M],£,(0)) = 0, 

and 



([M m ],[M],£,(0),£ z (0)) = -([M m ],£,(0),[M],£ z (0)) 
Finally, we see that 

(4-21) z m {s) = J] ( gmQ ^ ( ° )|2 (1 - 008(2*151,))) , 



l a 2 ml \xi(0)\ 2 . 



i=i 



m 2 



Integrating equations (|4.21|) . we get 



(4.22) 



Zm(s) 



^/ g m a^|£,(0)| 2 r S in(2s|fl|,) 



1=1 



m 2 v mi 

Lemma 4.5. Not all of horizontal curves are geodesies. 



m = 1,2,3. 



m = 1,2,3. 
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Proof. To prove this proposition we present an example. The curve 

(S^ CL\\S^ 
~2,s, —,s,0,... ,0, — — ,d,c 2 

is horizontal with ci,C2 constant. Indeed, 

*(•) =^ 5(M!«,i) = f (, 2 - £ + , 2 - £) = 
i 2 ( S ) =0, i(M 2 x,x) = ^(-s 2 - y + S 2 + y) = 0, 

ZaW =0, ^{M 3 x,x) = ^(-y + s + y- s )=0. 

From the other hand, the curve c(s) does not satisfy the system (|4.16j) . The system (|4.16j) gets 
the form 

1 = 2(an#i - a 3 i# 3 s - 0216*2) 

= 2(— an^is - C121O2S + a 3 i# 3 ) 

1 = 2(a 3 i^ 3 s + a 2 i#2 + «u^i) 
= 2(a 2 i^ 2 s - 0316*3 - anOis) 

for the curve c(s). Summing up the first and the third equation, and then, the second and the 
forth ones, we write the latter system as follows 

(2 = 4an0i 

= -40116*15 

1 = 2(a 3 i6> 3 s + 0216*2 + andt) 
= 2(0216*25 - a 3 i6* 3 - aiiflis). 

We see that the first and the second equations contradict each other. □ 

Lemma 4.6. A curve c is a geodesic for the group Q n if and only if 

(i) c(s) is a horizontal curve and 

(ii) c(s) satisfies c(s) = 2Mc(s) ; I = 1, . . . ,n. 

Proof. If a curve is geodesic, then it is horizontal by Lemma 14.31 Proposition 13.41 implies that 
the vector c is also horizontal: c = Ya=i^i-^i- Since x(s) = 2Mx by (|4.16|) . we obtain the 
necessary result. 

Let the curve c(s) satisfy (i) and (ii) of Lemma 14.61 The horizontality condition (i) of 
Lemma 14.61 can be written in the form 

1 / \ dH 

(4.23) z m = -[M m x,xj =^\x\ 2 Am + (M m x,£) = — , m = 1,2,3. 

as in (|4.13j) . We see that c(s) satisfies the equations of the second line of (|4.3j) . The condition (ii) 
of Lemma f4.6l admits the form x(s) = 2Mx(s) in the coordinate functions. Define the following 
curve 7(5) = (x(s), z(s), £(s), 9) in the cotangent space, where 

(4.24) ( = ^-iMi(s) with 6* = (9 1: 9 2 , 9 3 ) constant. 

The relations ()4.24j) imply the equations of the first and the last lines of (|4,3|) . Differentiat- 
ing (|4.24j) . we get 



x 1,,. Mi . s 1 2 



£ = Mi = Mx = -M(2£ + Mx) = M£ - -@ z x, 

2 2 2 2 2 
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by the condition (ii) of Lemma 14.61 Q4.24JI . and (|4,1U|) , Thus, j(s) satisfies the Hamilton 
system (|4.3j) . Then, the projection onto the (x, z)-space, that coincides with c(s), is a geodesic. 

□ 



5. Connectivity by geodesics. 

Let us ask in the following question. Is it possible to join arbitrary two points of Q n by a 
horizontal curve? A theorem by Chow ^SI gives an affirmative answer. We present a direct 
proof and calculate the number of geodesics connecting the origin with different points. We 
need the next simple observation. 

Proposition 5.1. The kinetic energies £ = ^\x\ 2 , £ m = ^\i\\ m are preserved along geodesics. 
Proof. In fact, 

d£ 

= (^A m x,A m x) = 2(A m x, M^4 m x) =0 
by Lemma 14.61 and property (|4.8|) of the matrices A4 m . The same is for £. □ 



5.2. Connectivity between (0,0) and (x,0), i/O. 

Theorem 5.1. A smooth curve c(s) is horizontal with constant z-coordinates Z\,Z2,z% if and 
only if c(s) = (aus, . . . , a^nS, z\, Z2, z%) with au E M and Ya=1 ^2k=i a \i ^ ®- ^ n °^ er words, 
there is only one geodesic joining the origin with a point (x,0). 

Proof. Let c(s) be a horizontal curve with constant z-coordinates z\, Z2, z%. Then z m = 
and (|4.21[) implies 

n 2 I • (r\\\2 

= z m = e m ^ aml [^ sin 2 ( S |^), m = 1,2,3. 
2 \ 6 \i 

We define by continuity 4p ^ = s 2 at \9\i = 0. Since the sum Yl?=i "^p^ s ™ 2 H^lO * s 

not identically zero we deduce, that 9 m = for all m = 1,2,3. The Hamiltonian system (|4,3|) 
is reduced to the next one 

' x =2£ 
=(M m x,£), m = 1,2,3 

£ =o 

k 0m =0. 

We see that £ is a constant vector. Taking into account that x(0) = 0, we get x(s) = 
(a\±s, . . . , ct4„s) with ay = 2^1- This proves the statement. 

Now, let us assume that c(s) = (aus, . . . , a^nS, Z\, z^, z^) with constant z-components. Set 
as = (aus, . . . , Q!4 n s). Recall, that (M. m a, a) = for any vector a = (an, ... ,04. n ) an d 
m = 1,2,3. Then, 

1 • s 

i m = = -(M m (as), (as)) = -(M m a, a), m = 1, 2, 3. 

The horizontal condition (|3.1|) holds for all three z-components. □ 
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5.3. Connectivity between (0,0) and (0,z), z / 0. We need to solve the equation (|4.16f) 
with the boundary conditions 

x(0) = x(l) = z(0) = 0, z(l) = z. 

We also need to know the initial velocity x(0) since we do not have enough information about 
the behavior of x-coordinates. In the following theorem we use the notations n = (m, . . . , n n ), 
rii G N, for I = 1, . . . , n, and 

r l 



(5.1) 



N 



o 







-U 



Theorem 5.2. There are infinitely many geodesies joining the origin with a point (0, z). The 
corresponding equations for each n = (m, . . . , n n ) are 



(5.2) x{ n) (s) 



1 — cos(2s7rnz) 



[Z]jii(0) + 



sin(2s7rn;) 



(irni) 2 L 11 iy ' ' 2-kui 
where Z is a block diagonal matrix with the blocks 



U±i(0), 1 = 1, 



(5.3) 

and 
(5.4) 








23«3I 


22«2i 




l*(0)l^3 


I*(°)ISm 2 









23(13; 


\m\* NAl 


\m\%A 2 


l*(0)lk 3 




22<l2i 





21 an 


l*(0)l^ 3 


l*(0)I^A 2 




22^2! 









l*(°)l^ a 







^ ( s ) 



|i(0)|f 



E 



a m«l X 



IJV4 m i = i 

The lengths of corresponding geodesies are 

3 4(1) 
m=l 2^/=l 



sin(2s7rnj) \ 
27rn; / ' 



m = 1,2,3. 



; 2 = 16 V ml J 



16 E 



zL(l) 



Proof. Substituting s = 1 in ()4.20|) . we calculate 



i 1^(0)1 



= |x(l)| 2 = J>(1)| 2 = £ 



z=l 



sm z \6\i 



E 



xK0)| 2 . 



1=1 



Since the kinetic energy £ = l^ML does not vanish, there are indexes I, such that |x;(0)| 2 7^ 0. 
We deduce in this case that 



(cos2|0| 



A\e\f 



[M],i,(0), [M]/xj(0)) + 



sin 2 2|6>b , . 



|0| z = y 0f a\ + 0|ct| + 0|a| = nn u n x e N. 
If |x/(0)|=0, then the corresponding \6\i are arbitrary. Equalities 1)4. 22j) give for s = 1 



(5.5) 



z m (D 



a^|x,(0)| 2 _ m ,^ |2 



^ 4(vrn;) 2 



|x(0)| 



NA m 1 



m = 1,2,3. 
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We find the unknown constants 6 m = i-t^T^ ■ Substituting 6 m in (|4.2U|) . (|4.22|) . we obtain 
the equations (|5.2j) and (|5.4|) for geodesies. 

To calculate the length of geodesies, we observe that Ylm=i @m a mi = 7T ' 2 " n 'i an d deduce 

^ z m {l)B m - - ^ ^ T^vj - 4 2. T^-p ^ a ™< - 4FWI - 2 

m=l m=l Z=l v ; Z=l v ly m=l 

from (|5.5|) . The lengths of geodesies are 

^ > - 4(i) 



ll = (J q \x(s)\d s y = \±(0)\ 2 = A Y^z m 9 m = 16 



, v n a ^l A '(°)l 2 

m=l m=l 2^2=1 (Trn;) 2 



Z 2 (T) 

(5.6) = 16 ^ " 2 7 , n = (m,...,n n ), n, e N. 

m=1 Fl u JljVA m 



□ 



Remark 5.3. Let us discuss the cardinality of the set of geodesies. In the general case in Theo- 
rem E21 when a m i are different we obtain the countably many geodesies connecting the origin 
with a point P = (0, z). If the multiindex n = (m, . . . , n n ) increases, then the geodesic rotates 
more frequently around straight line, connecting origin with P = (0, z) approaching to the line 
and in the limit we obtain a limit curve of Hausdorff dimension 2. We present a graphic of three 
geodesic curves (|5.2[) . (|5.4[1 with the initial velocity x(0) = (in, 0, 0, 0, £21, 0, 0, 0, 0, . . . , 0) and 
the end point P = (0, 21, 0, 0), where in, ±21,21 are different from zero. The corresponding 
multiindexes are n = 1, n = 2 and n = 5. 




Figure 1. The graphs of geodesies where the vertical axis represents z coor- 
dinate and horizontal exes are x\\ and X21 

The projection of geodesies into the horizontal subspace belongs to an ellipsoid passing 
through the origin. 

If, in particulary, an = 021 = 03; = o; then rearranging the indexes, we can assume that 
«i < 0-2 < ■ ■ ■ < cLp = Op+i = ■ ■ ■ = a n . Applying the rotation to the geodesies in the subspace 
(0, . . . , 0, x p , x p+ i, . . . , x n , 0, 0, 0), we get uncountably many geodesies. In this case we have the 
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following estimate of their lengths: 

/ 



2 _ 16|*(1) 



1^1=1 (7rn ; ) 2 



If a m \ = a m 2 = . . . = a mn = a m , then the multiindex n reduces to the index k £ N. The 
equation (|5.6jl implies 



^ = |£(o)|* = 4vrA: 2J J! y i , A; G N. 

m=l 



Let £7 be a neighborhood of the origin O. From Theorem 15, 21 we know that no matter how 
small U is, we can always find points in U which are connected to O by an infinite number of 
geodesies. This is totally different from the Riemannian geometry. It is known that every point 
of a Riemannian manifold is connected to every other point in a sufficiently small neighborhood 
by one single, unique geodesic. 

5.4. Connectivity between (0, 0) and (x,z), x ^ 0, z ^ 0. Now, we will look for a solution 
of the equation (|4.16|) with the boundary conditions 

x(0) = 0, z(0) = 0, x(l)=x, z(l) = z. 

Let us make some preliminary calculations. We obtain 

(5.7) |±K0)| 2 = _i|L|^(i)|2 ; / = i,..., n , 

sin 

from (|4.20|) for s = 1. Putting s = 1 in (|4.22|) and making use of (|5.7j) we obtain 

(5-8) ^(l) = %E ^ l |g | (1)|a M(|g|0, - = 1,2,3. 

i=i 1 

where fi(\0\i) = . , , — cot(|0|z). The function /j.(9), introduced by Gaveau in [T^J, was first 
studied in detailed by Beals, Gaveau, Greiner in O [7] - By the following lemma, one finds 
some basic properties of the function \x. 

Lemma 5.5. The function fi(9) = ^ e ~ c °t is an increasing diffeomorphism of the interval 
(— 7T, tt) onto R. On each interval (run, (m + l)vr), m = 1, 2, . . ., the function \x has a unique 
critical point c m . On this interval the function \i strictly decreases from +oo to /i(c m ), and 
then, strictly increases from n(c m ) to +oo. Moreover, 

M(cto) + vr < /x(c m+ i), m = l,2, 

The graph of fj,(9) is given in Figured 

Theorem 5.4. Given a point P(x, z) with x\ ^ 0, I = 1, . . . , n, z / 0, there are finitely many 
geodesies joining the point O(0,0) with a point P. Let = (|#i|i, • • • , |#i|n)j • • • > $(N) = 
(\6n\i, . . . , |0Ar| n ) be solutions of the system 
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Figure 2. The graph of n{9) and solutions of the equation /u(|6>|) = 



VFe fix one of the solution •& = (|0|i, . . . , \0\ n ). Then the equations of the geodesies are 
(5.10) 

x[ n) ( S ) =(4cot(|0|Osin 2 ( S |0|O-2sin(2 S |0| / ))|^x / (l) 

+ (- cot \9\i sin(2s|<9|0 + sin 2 (s|0| i ))^x / (l), l = l,...,n, n = 1, 2, . . . , N, 



n "LMWfi f„ gin(g£lglj 



(5.11) z£>(a) = zm(l) 



En 
1=1 



)) 



, 771=1,2,3, n = 1,2,..., N, 



where Z is a block diagonal matrix with blocks (j5,16|) . The lengths of these geodesies are 

4(i) 



(5.i2) , l= ±m>-w 



i=i 



16 Y, 



Proof. We have 
(5.13) 

from (|5.H|) . Then 



En 
Z=l 

4z m (l) 



K\o\i) 



+ X>(i)| a |0|,cat(|0|,). 



1=1 



, m = l,2,3, 



^ z m a ml 



3 3 

(5.14) |0| 2 = y 9la 2 ml = y —m-rru Z = 1, . . . , n, 

^ ^i ^ WIPA) ^ 2 



m=l 



that prove (|5.9|) . 

Let us fix one of the solutions of the equation Q5.9|) i? = (|0|i, • • • , \0\ n ) for a given point 
P(x,z). Putting (f5~7)) and (jSHSj) into (g3Z|> , we obtain (|K7TT|) . 
Setting s = 1 in we find x[ n) (0) for = (|6>|i, . . . , |0| n ): 

x[ n) (0) = 2|0| z [sin(2|0| z )^+ (l-cos(2|0| i ))^]~ 1 xKl) = [(\e\i cot - [M] J 
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This and (|4~2U)) imply 



(5.15) 



= I[(2cot|e| z sin 2 ( S |0| / )- S in(2 S |0|O)^ 
+ ( cot \6\i 5111(25^1/) + 2 sm 2 (s\9\i))u} xi(l), 



1 = 1, 



, n. 



Taking into account (|5.13|) . we deduce that each block [M]/ of the matrix M takes the form 
4[Z]j with block [Z]j written as 



(5.16) 







~2T — i2~~ 



Ef =1 ^f^Mieh) Ef =1 ^f^MI^I;) ElLj ^^(l«IO 



23(1)03; 





22(1)02; 



22(1)02; 



23(1)03; 



ELa ^f^(|e|;) Ef=a ^^Ml^l;) 



|0|; 
22(1)02; 



|S|; 
23(1)03; 





21(1)01; 



2i(l)oi; 
2 i T .|2 



Ef=i ^f^Miei;) 




191; 



Finally, (|5TT5|l and (f5~T5|) give (|5~TU|) . 

To obtain the length of the geodesies, we make the following calculations. 



3 1 n 



m=l 



Z=l 



N(i)M^) 



3 1 n 



771=1 



(5.17) 



l^ |xKl)l 2 |g|f i r , m|2|fll ., lfln 
li(O)] 1 



sim^l/J 



i=l 



1 n 

-^|x/| 2 |0bcot(|^). 



From the other hand (|5. 13f) implies 



(5.18) 



E *»d)*». = 4 E — 



m=l 



-1 E^^-S^Mdei,)' 



The formula (|5.12|) follows from (|5.17|) and (|5.18|) . 

Remark 5.5. Let us consider the particular case an = 02/ = 03; = a/ > 0. We have 



□ 



(5.19) 



4\z\=J2MMi)\ 2 »We\) 



1=1 



from (|5~H1) . Here \6\ 2 = 61 + 6% + 6%. We denote by |0|i, . . . , \6\ N the solutions of (OH) and 
let |0| one of the solutions. Then (j5,13|) implies 

4z m (l)|fl| _£m |fl| 
E?=ia*M 2 MN0|) N 
To obtain the formula for the length of geodesies, we write (|5. 19|) as 



1 - 



sin 2 (|0L) 



^|a i |N| 2 cot(|0| / ) = ^-^|a i |^| 2 cot(|0| / ) 



z=i 



Z=l 
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and get 

n 

ll = \9\{A\z{l)\ +^a i |x / (l)| 2 cot(a^|)). 
1=1 

Simplifying more the situation and supposing that a; = a > for all / = 1, . . . , n, we get 
that \6\i = a\9\. This implies that \9\ is a solution of the equation (see Figure |5J) 

(5.20) »(a\9\) = ij^L 

In this case to calculate the length of geodesies joining (0, 0) to (x, z), x\ ^ 0, I = 1, . . . , n, we 
use the homogeneous norm \(x, z)\ 2 = \x\ 2 + i\z\. It gives for a solution \9\ a , a = 1,...,N, 

of EH 

| x |2 + 4 | z | = | x |2 + a \ x \i^ a \e\ a ) = (1 + Q j »(a|g| a )) Sm 2 ( " |g 2 la) /2 
by ((5371 and ([5~7fl . Then 



2 ° \9\ 2 

a sin(a|#| Q ) ( sin(a|#| Q ) — cos(a|#| a )) + a 2 |#| Q ^ +11 

In the last simplest case it is easy to observe that if z is fixed, and |x| tends to zero, then the 
ratio increases and the number of solutions of the equation 4^ = fj,(a\9\) also increases 
(see Figure EJ. In this case, the function fj,(a\9\ a ) = a|e| " cos MeU)sin(a|e| Q ) tendg tQ infinity ag 

sin (a|p| a J 

|x| — > 0, and we obtain that sin 2 (a|#|) = and a\9\ = nn, n £ N. One sees that Theorem 15.21 
is the limiting case of Theorem 15.41 as the ratio tends to oo. If we fix x and let \z\ tend to 

0, then the equation 4^ = fi(a\8\) says that n(a\9\) —* 0. This implies that \9\ — > and we 
obtain Theorem I5.1l as another limit case of Theorem 15.41 

The last particular case is when a m \ = a m i = . . . = a mn = a m . We denote \9\f = a\9\ + 
a^9\ + a\Q\ = \9\ 2 a . The equation to find \9\ a is 



K\9\a) 



The value of m and the lengths of geodesies are 




*m(l)\e\* l 2 (\6\ a ) = \9\j4 i Z M± + \x(l)\ 2 cot( 



,2 



^m(l) \ V O: 



In the following theorem we consider the connection between the origin and a point P(x, z) 
when some of the coordinates x; vanish. 

Theorem 5.6. Given a point P(x, z) with xi ^ 0, I = 1, . . . ,p — 1, and X[ = 0, / = p, . . . ,n, 

z 0, there are infinitely many geodesies joining the point O(0,0) with a point P. Let S\ m = 

XSli aroj |gj t (1)l K\9\l), S 2m = Ya= p ^t^S^ ' n l 3 = ( n P' ■■■■> n n) be a multiindex with positive 
integer-valued components for each (5 € N, and "& K = (|0|i, • • • , k = 1, . . . ,N be solutions 

of the system 



(5.21) WS-± ^1% , ( = 1,. 
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Then the equations of geodesies are 
(5.22) 



x[ K \s) =(4cot(|0|z) S in 2 ( S |0| / ) -2 B m(2 a |0|,))^x,(l) 



- cot \6\i s^l^lO + sm 2 (s\9\i))Uxi(l), l = l,...,n, k = 1,2,...,N, 



(5.23) 

where 

(5.24) 



x^\s) 



( 1 — cos(2s7rn/) 
(vrn;) 2 



Z]jxK0) + 



sin(2s7rn;) 
2-kui 



U±i(0), l = l,...,n, /3eN, 





S11+S21 
S13+S23 

S12+S22 



S11+S21 



Z2Q2i 
S'l2+S'22 

S13+&3 



2303; 

5'l3+523 
2202I 



■Z2Q2i 

S12+S22 
Z3«3i 



S12+S22 S13+S23 

Sii+S^i 
21Qli Q 



511+521 



and 
(5.25) 

(K,n^) 



(1) /v ^(l)l t V gnggjgji) . 



fo(o)l 5 



Sim + ^ 7^ Sin 



7r 2 n? 



sin(2s7rn/) 
27rn; 



wrf/i m = 1, 2, 3. 

T/ie lengths of these geodesies are 



(5.26) 



/ 2 = V 



Xl 



1=1 



16 E 

771=1 



p-i 



Sim + S 1 ; 



2f» 



+ X)l^(l)| a |«|ioot(|fl|,). 



Z=l 



Proof. If xz(l) = 0, then the formula 



sin 2 (|e| ; ) 



I ±1 (0)1 2 implies that \xi(0)\ = or 



sin 2 (|#|;) = 0. If |i;(0)| = 0, then the corresponding xi(s) = 0. The more interesting case when 
\±i(0)\ ^ for I = p, . . . , n. Then \9\i = irui, rii € N, I = p . . . , n. We deduce from (|4~2^|) for 
s = 1 

(5 . 27) . ^-^(^aUxm?,,^ , ^c&|s,(0)|*' 



, m(1) ^ (E «^ ) + E , 

;=i 1 " z= P i 



— {Sim + ^m)) 



where the number ni can have any positive integer value. We conclude that the sum S\ 



EI 



" "ml 



<«l*l(0)l 2 



admits countably many values. To define we find 6 n 



4z m (l) 

5l m +S2 ri 



2 in 



m 



1, 2, 3, from (|5.27f) and then argue as in Theorem 15.41 

3 3 
(5.28) |#| 2 = £ 0^ = E 



!6^(l)aL 



i = l,...,p-l. 



777=1 777=1 (5*1777 + ^m) 

Conclude, that for each multiindex with positive integer- valued components = (n p , . . . , n n ), 
j3 G N, the equation (|5.28|) defines the multiindex = (|0|i, . . . , |0| p _i), k = 1, . . . , N. Let us 
fix one of the solutions ($ K ,np) = (\6\i, ■ ■ ■ , \0\ p -i,n p , . . . ,n n ). The relations (|4.2()|) for s = 1 
give 

i,(0) = fdOljcot - pM0,V*(l), ? = 1,...,P-1. 



(5.29) 

Substituting (t5~29l) into (lOTl and (E32l) . we obtain (t5~22l) . (Q3|l . and (t5~25l . We get the 
form of (|5.24|) from 6 



iz m (l) 



m - Si m +S 2m ' 



m = 1,2,3 and the definition of the matrix M. 
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To calculate the length of the geodesic we argue as follows: 



' ' W n 2 V^I^(1)| 2 MM/) 
/ j "rn a ml / j 



m=l 



m=l 

p-1 



1=1 



\0\i 



Xl 



2^,2 



m=l 
2 



i=l 

p-1 



i^ \ Xl (i)\ 2 \e\ 2 i p ' » |^ (0 )p 

iL-^T-iLW 1 )! l*li«*(l*li) + 2^— r" 



(5.30) 



4^ sin 2 (|^ 

|2 i P' 1 

~ 4 

«=i 

4z m (l) 



1=1 



l=p 



■r 



1 

-J2\^)\ 2 \0\icot(\9\l) 



From the other hand, since 9 m = -s — "r Q , m = 1, 2, 3, we deduce 



Sim 4" 02m 



(5.31) z m m = 4 ^ 

m=l m=l 

The formula (ETSfill follows from (t5~3Ul) and (153T1) . 

Remark 5.7. Let make some simulations for the anisotropic group Q 2 . Set 

xi(l) = (x n ,x 12 , 0, 0), x 2 (l) = 0, x 2 (0) = (i 2 i(0), i 22 (0), 0, 0), 
zi(l)^0, z 2 (l) = z 3 (l) = 0. 
In this case the equation (|5.21|) can be written in the form 



(5.32) 



Mlfli) 



4Ui| 



L af 2 |x 2 (0)| ; 



□ 



an|xi(l)| 2 7r 2 n 2 af 1 |xi(l)| 2 ' 

We present the solutions for different values of n: n = 1,2,50 in Figure El We see that for 
sufficiently big value of n the second term in the right hand side of (|5,32|) goes to 0, and we 
obtain a finite number of solutions for 101 1. 




2.5 5 7.5 10 12.5 15 17.5 



Figure 3. Solutions of equation (|5.32|) 

Moreover, we obtain countably many geodesies, because of the second part of multiindex, 
corresponding to the positive integer values is countably infinite. Nevertheless, since the sums 
S2m tends to and the sums S'lm are strictly positive as n — > oo, we conclude that the lengths 
of these geodesies are bounded from the above. The projection in each subspace x\ are still 
ellipsoids. In Figures 0] and El we present the projection of a geodesic into spaces (xi,zi) 
and (x2,zi). We can see that the number of loops is different and increases in the subspace 
corresponding vanishing value of xi(l). 
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Figure 4. Projection of a geodesic to the space (xi,z\) 




FIGURE 5. Projection of a geodesic to the space (2:2, z\) 

6. Complex Hamiltonian mechanics 

Our aim now is to study the complex action which may be used to obtain the length of real 
geodesies. 

Definition 6.1. A complex geodesic is the projection of a solution of the Hamiltonian sys- 
tem Q4.HJI with the non-standard boundary conditions 

x(0) = 0, x{l)=x, z{0)=0, z(l) = z, and 

8 m = iT m , Tfl = 1,2, 3, 

on the (x, z)-space. 

Let us introduce the notation —it for the vector (— ir%, — IT2, — IT3). We write |r|/ = 

\j a li T l + a l r 2 2 + 4 T l- Then \°\i = \] a li 9 l + a li e l + a li e l = i\ T \i- 

Notice, that we should treat the missing directions apart from the directions in the under- 
lying space. 
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Definition 6.2. The modifying complex action is denned as 
(6.1) f(x, z, t) = -i V" T m z m + / ((&,£) -H(x,z,£,t)) ds. 

We present some useful calculations following from the system (j4.3j) . 

= m 2 + (Mx,o = i|x| 2 - ^(Mx, x ), 

(6-2) |£| 2 =M! _ I(Mx, x) + J(Mx, Mi) = ^- i(Mx, x) + ^(0 2 x, x), 

(Mx,0 =^(Mx,±)-^(@ 2 x,x). 
Making use of the formulas (|4.2|) , (jo,2j) , and (|5.7|) , we deduce 

f(x, z,r)=-tV r m z m + / ((x, £) - (x, z, £, r)) 

= -i^T m z m + I (MfL - i(Mx,x)) ds 

m J° 

l^(0)| 2 f 1 , , , . N , 
'2^ r m 2; m + 2^ — 4 — y cosh(2s|r|/)ds 



I 

m 1=1 

E^. bzl 2 . , . 

TmZ m + 2^ — 2~" l T U coth \ T \l- 



V- , ^ \xi\ 2 {i\r\i? sinh(2|r| ; ) 

2^r m z m + 2 ^ 4 ^(.j^i,) 2 |r|, 



4 



The complex action function satisfies the Hamilton-Jacobi equation 
(6.3) £ r m |^ + H(x, z, V x f, VJ) = f. 



Indeed, we have 



m=l 



^p- = -iz m - %T m arr f) X } fi(i\T\i), m = 1,2,3. 



lx,z '3x'^ j a ^ z ^ T > 2-, 4 sinh 2, r 



from g2J), (fo~2j) . and g^SJ). Then, 



3 o /• ^ <■ o ^ n , , 2 | 



m=l m ;=1 



^ + ^2 ~^-\ T \l C ° th l T l« = /• 
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In the critical points r c , where -qJ- = we have from 1)5.1 7jl 

f(x, z, r c ) = H(x, z, V x f, V z f) = | = ^( 7 ), 
where a geodesic curve 7 connects the origin with (x,z). 

7. Green's function for the Schrodinger operator 
Consider the Schrodinger operator 

r a ■ 9 

ou 

We are looking for a distribution P = P(x,z,u) on x M.^, x M+ satisfying the following 
conditions 

1) LP = A P - i§£ = for u > 0, 

2) lim P(x, z, n) = 5(x)(5(z), 

u->0+ 

where 5 stands for the Dirac distribution. 
The next propositions are easily verified. 

Proposition 7.1. For any smooth function <p and any smooth vector fields X\, . . . ,X n we 

have 

Ae* = e' p (Aip + \V<p\ 2 ), 
where A = E]=i X j and l Vl ^l 2 = £?=i(*j¥>) 2 - 

We recall that X denotes the horizontal gradient Xn, . . . , X± n . 
Proposition 7.2. 

4 n 

H(x, z, V x f, VJ) = \Xf\ 2 =Y,J2 ( X uff- 

fc=i 1=1 

Proposition 7.3. Let V and f be smooth functions of x, z, t, and X±, . . . ,X n smooth vector 
fields. Then for any number p, we have the following identity: 

(7.1) A(vrn = (Avor p - vr p - 1 [W)v + 2(v/)(w)] + (- P )(- P - i)/-*>-viv/i 2 . 

Proof. The formula 1)7.1)1 is obtained by the direct calculation. □ 

Before we go further, let us make some calculations. We apply Proposition 17.11 to e « and 
system 1)2.3)1 of horizontal vector fields X^, fc = 1, . . . , 4, I = 1, . . . , n. Introducing the notation 
tp = -if, we get A </? = -£A /, \Xip\ 2 = -^\Xf\ 2 , and 

* • - 1 _ ^V(r)u 2n ^( ; . . 1 . |2 



(7.2) A e*> = e*>l- -A f - -\Xf\ 2 ) = ^' - zA f - -\Xf\ 

V u u / u^ n+ V{t) V u 

Hamilton- Jacobi equation ()6.3)) . Proposition 17.21 and the equality Q7.2JI imply 
( 7 - 3 ^ A o^r = { ~ ^0/ " - + - 

m=l 



Differentiating e J^ff with respect to u, we obtain 

9 fe^V(r)\ e^Vir) (f . ,\ 
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Summing (|7.3|) and 1)7.4(1 . we have 

m=l 

We express — ^ ^m=i T ™-Th~ from the formula 

E ^ (e»V(r)r ra ) = e-V(r) ( - 1 £ r ra ^-) + + 3e^(r) 

m=i aTm m=l aTm m=l ° m 

and put it into (|7.5|) . Finally, we deduce 
(7.6) 



m=l 

3 



+^E^(^M 



u 2n+4 Q r 
m=l 

The equation 

3 BV 

(7.7) (2n - Af)V - E ^m*— = 

m=l 

is called the transport equation. We show that the function 

n I |2 

^)=Il7wr 

7-7 sinn r ; 

i=i 1 11 

is a solution of transport equation. Indeed, since 

" N 2 



we have 



/ = f(x,z,r) = -i^r m z m + ^ ^p|r|/coth(|r| 

m 1=1 

df d 2 f 

-iT m , TrT - = 0, m = l,2,3. 



o "mi ^2 
OZ m UZjn 

1T~ = \x k i\T\icoth{\T\i) } = ^-c5oth(|r|j), fe = 1, . . . , 4, Z = 1, . . . , n. 

ox k i 2 dxfa 2 



Finally, 



and 



A/ = 2^|T| iC oth(|r| 



l=i 



(7.8) (2n - A/) y(r) = 2V(r) (n - ^ |r|, coth(|r|,)) • 

i=i 

On the other hand the equalities 



dV A A \r\-f d ( It 



e n 



9r m ^ z= f^ r sinh 2 (|r| z ) 9r m V s inh 2 (|r 



2 



Irl 2 2« 2 r 



EII S^fe-S^d-lrtW)). n.- 1,2,3, 

r=ll=l,l^r Vl |ly 1 
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imply 

m=l r=li=l,t^r vl 1 ' vl 1 ' 

n I |2 n n 

= 2 II ■ lln i \ • " M' coth (Mr)) = W(T)(n - \r\r coth(|r| r )), 

1=1 SmJl vl r Ui r= i r= i 

that shows that V(r) is a solution of the transport equation ()7.7j) . The function V(t) is called 
the volume element. 

If the volume element V(r) satisfies the equation ()7.7|) then the equation (|7,6|) is reduced to 
the next one 

m=l 

We note that the expression e v V{T)r m vanishes as |r| — > oo. Integrating over R 3 with respect 
to dr = dTidT2dT3, we obtain 



( A ° 1 du) J m 



dr = for u > 0. 



Thus, the function 



c r =it 

P(x,z,u) = T / e u F(r)dr 

satisfies the first condition to the Green function at the origin of the Schrodinger operator. 

7.4. The heat kernel on Q n . In this section we denote the time variable by t and we will 
consider the heat operator 

A o-|=E^-|- 

k,l 

where Y = (Yn, . . . ,F 4n ) = V y + ±(Ylm=i M ™ys§^) with y = ( yil ,.. . ,y 4n ). The fundamen- 
tal solution at the origin is the function P(y, w, t) defined on Q n x such that the following 
conditions 

1) A P - = for t > 0, 

2) lim P(y,w,t) = 6(y)8(w) 

hold. With the change of variables 

it = it, x k i = iy k i, w m = z m , m= 1,2,3, fe = l, ...,4, Z = 1, ...,n, 
the heat operator transforms to the Schrodinger operator 

du ^ kl du 

k,l 

Indeed, under this change of variables we obtain 

d d 

— i— = — and X = —%Y. 
ou at 



The calculus of the previous subsection give us the following statement. 
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Theorem 7.1. TTie /leaf kernel at the origin is given by 



C f -f 



where 

n I |2 

f{y,w,r) = -i^TmWm + ^ |r|j coth(|r|i) 

m 1=1 

is the modified complex action and 



II 

v{t) = n sinh 2 ( | T | 



zs i/ie volume element. 



7.5. Green function for the sub-Laplace operator. Let us integrate the kernel P(x, z, u) 
with respect to the time variable u on (0, oo). That is 



P(x, z, u) du 



POO (~t r 

~- / -T-^ / e- tS ' u V(T) dr du 

--C [ V(t)( f°° u- 2n - 3 e- if / u du)d,T. 
Jm.3 Wo 1 



We first look at the inner integral: 



oo 

u -2n-3 e -if/u d ^ 



10 

Changing variable v = ^, u = ^, yields dv = = Zf du and du = dv. Hence, we have 



u 

oo i roo 



2n-3 -if/u d = i / -v 2n+3-2 rf = T(2n + 2) 

C j2n+2j2n+2 J Q j2n+2 pn+2' 



2 f 2n + 2 l 

Let us introduce the following notation 

~G(x, z)=J q P(x, z, u) du = C T -^±fi p Z { ^ ZjT) 

The aim of this section is to show that the function —G(x, z) is the Green function for the 
sub-Laplacian operator. Firstly, we need some auxiliary results. 

Proposition 7.6. Denote 

n -y 3 3 

f(x,z,w) = ^2 -|x/| 2 H/coth(H/) - i ^2 w m z m = l{x,w) - i ^ w m z m , 

1=1 m=l m=l 

where w = r + isz. Then there exist positive constants c\, c 2 , and eo such that for all real T m , 
all < e < eo, and all x 6 M 4n ; z = (z\, z 2 , £3) € R 3 we have the estimates 

(7.10) \lm(i)(x,T + i£z)\ < cie\x\ , 

(7.11) Re(j)(x,T + iez) > c 2 \x\ 2 , 

(7.12) Re(f)(x, z,t + isz) > c 2 (\x\ 2 + e\z\). 
Here £ = jfj */ z 7^ anc ^ z = if z = 0. 
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Proof. If z = 0, then Im( 7 )(a;,T) = and since |r|/ coth(|r|/) > 1, I = 1, . . . , n, we have 
Re( 7 )(x,r)>M!. 

X^ m= i a ml( T m + i£z m ) J = ai + iPi, where 

\ i /a circtcin ( — — "9^1 ~ 1 2 ~ — ) 
r|?-e 2 |5|2) 2 +( 2e ^aLr m 5 m ) 2 ) cos V 7 + nd, d = 0,l, 



ot-l = 

m 

and 



. 1/4 arctan , , -. 

A= ((|r| 2 -e 2 |z| 2 ) 2 +(2£^aLr m 5 m ) 2 ) sin V f f 7 + vrd, d = 0,l. 



We consider the case d = 0, another one can be treated similarly. Since 

. . n . sinh 2a sin 2/3 

coth(a + ip) = — — — — 1- 



cosh 2a — cos 2/3 cosh 2a — cos 2/3 
we have 

™ Ix I 2 

Re( 7 )(ar, r + iez) + ilm( 7 )(x, r + iez) = ^ — ^— (a/ + ify) coth(ai + ifa) 

1=1 

I xi 1 2 / a\ sinh a; cosh a; + (5\ sin /?/ cos ft, 



E 



, . 4 V ginh^ ai + sin 2 /3/ 

I x/ 1 2 / Pi sinh a; cosh a/ — a; sin /3/ cos Pi 



■ sr^ \ x i\ ( Pi smn a i cosh ai — ai sin p; cos p/ \ 
1 4 V sinh 2 a^ + sin 2 A / ' 

Denotes by ipi the angle between nonzero vectors (auTi, (I21T2, (I31T3) and (auzi, CL21Z2, CL31Z3). 
We consider two cases, when cosipi = for all I = 1, . . . , n, and cosV'z = 1?/ ^ for some 
index Z. 

Case 1. If cos ipi = for all I = 1, . . . , n, then ^™=i a mi T mZ m = 0. We have 

a,= (|r| 2 - £ 2 |5| 2 ) 1/2 and /?, = 0. 

It gives 

A> |x/| 2 / A sinh q ; cosh a/ - a; sin /3/ cos /3/ \ 

Im( 7 )(x,r + iez) = > — — -= =0; 

^ 4 V S inh 2 a ; + sin 2 t ' 

Re( 7 )(ar, r + iez) = —^-ai coth(a;) > — ^— , Vq; G E, 
z=i 

because a/coth(a/) > 1. 

Case 2. If cosipi = "&i 7^ for some / = 1, . . . , n, then Ylm=i a mi T mZm = $i\t\i\z\i. We can 

I 1 2 

suppose that e satisfies < e 2 < min; =lj . . „ {^p}- We obtain 



\r\, < HF^W Ml ' 

(7.13) fcxMK (^ + (2^|rM%) 2 ) 1/4 < ((M 2 " ^l^?) 2 + (2^E m «Lw m ) 2 ) 1/4 

< (|r|f +(2^|r| / |5| i ) 2 ) 1/4 <fc 2 ( 1 9)|r| i . 
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Now, we put one more restriction to s assuming that e < mim = i r .. jn { ijngn }■ Then 

2e^i\z\ t 1 2e0,|z|, 1 2ct?,|t|,|z|, . 1 . 4ei9 { |«|i . 2e0,|5|, 

— — : — < — arctan — : — : — < — arctan - — p= 9 9 < — arctan — |— : — < — : — : — , 

A T \i 2 \ T \i 2 \ r \i- £ i\ z \i 2 \ T \i \ T \i 

and we get 

arrtar, ( ^WlWl ) 



(7.14) ^ < cos j < cos ^ < cos j < 1 



arctan f 3s fll r lil £ lj ', 

< sin < sin -BlHEH^. < sin < H« 



\^F\rj ^ 0111 2 " lzli ' - — v^^y " "Fit 

We observe that |z| 2 = X^m=i a miT^Y I — Hl=i a mZ — ^> where a = max{a 2 ,}. Prom the other 

II m,i 

estimate the value of a; and $ as follows 



hand, if we denote a = min{a^}, then a < Yl m =i a mifffi = \^\f- F rom <|7.13|) and (|7.14|) we 



(7.15) h\r\i < ai < k 2 \r\i, 



(7.16) h(a)e < Pi < fc 4 (a)e 

If \t\i < 1 we use the Taylor decomposition and obtain 
Pi sinh ai cosh a; — ai sin Pi cos Pi 



sinh 2 a; + sin 2 Pi 
If |t|z > 1 we argue as follows 

Pi sinh ai cosh a; — a; sin /?/ cos Pi 



laiPi(af - Pf) + 0(af - Pf) 



af + pf - 0(af + Pf) 



< he. 



sinh 2 ai + sin 2 $ 



< 



coth(a;)| + 



because ai is bounded from below, the functions |coth(a/)| and 
above. The last two estimates imply 



sinh (cti) 

I Qi I 
I sinh 2 cti I 



|) < k 6 e, 

are bounded from 



1=1 i 12 

xr 



|Im(7)(x, r + iez)| < ^ —^k-je < c\e 



To obtain (|7.11|) we change the arguments. Let us focus on the value of the derivatives 
dl Qw^ at w m = Km, Cm G K for m = 1, 2, 3. The equality 



d"f(x, w) 



dw r , 



Ti i 1 2 2 

X U a ml W m 



w=iC, 



E 



w|z V s inh 2 (|w|;) 



coth(|u;|;) 



n I 1 2 2 a 
;5-jfk^ M( | C | {) 



Z=l 



Kb 



implies ^g£^ 
Let us look at the 

a 2 7 



= and we conclude that w = iQ is a critical point for Re(7(x,w)). 
dessian at w = iQ. We have 



E 

l=i 



\ X \f a ml 



Kb 



a 2 (^ 2 

_miym 

ICI 



2« 2 r 2 

2 , • ,,, 2 Pferr Ci-IChootacii)) 



) + 



lC| 2 sin 2 (|CI0 
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Since 1 



a 2 C 2 

"mi W 



> and 1 - |CU cot(|CU) > we see that 



a 2 7 

dw 2 



m = 1,2,3. The mixed second derivatives are 



> for + |C|/ < f , 

■w=iC, 



5 2 7 



dwm dwv 



M^I 2 «L«LC m arcot(Kb) 2|cUcot(KU) 



w=iQ 



E 

1=1 

n 

E 

1=1 



ICI? 



k/l 2 4/4CmCfc 



Id? 



ICU 



+ 



sin 2 (|C|0 sin 2 (|C|0 



-. We observe that since all second derivatives of 



where 5 (|C|/) - S^^T + itaW 

7(x, w) are real at the critical point, the Hessian for 7(3;, w) coincides with the Hessian H for 
Re(7(x, w)) at w = iC- We write M-Z- as 



5 2 7 



dw"L 



w=iQ 



E^ 



z=i 



4 L 



ICI/ 



+ 



ICI? 



l -g(\(\i) 



Then the Hessian can be written in the form H = Ya=i Hi, where 



Hi 



"?,MICI0 , af,C : 2 



^fpsGCIz) 



a li a liClC3 



icTT 



*(ICIi) 



To show that is positive definite we need to show that each Hi is positive definite. It was 
shown that 



L ^(ICU) 



4MICI*) , 4Ci 2 



+ 



Then we have 

4MICI/) . 4d 2 
ICI 2 



ICI/ ' ICI 2 

4m(ICIz) , 4cf 
ICI 2 

,2 „2 



3 (|C|/)>0, for 0^|fl/< 



7T 



<KlC|l))( 



+ 



g(\C\i) 



ICI 2 



ICI/ ICI 2 / V ICU 

44m 2 (IC|/) (, 4C 2 + 4Cf \ , 244MlCl/)(4C 2 + 4C 2 2 ) h ... rn+nmW n 
^ i 1 ^ J + |C|fsin 2 (|C|/) (1 IC|iC ° t(ICI ^ > ° 



ici 2 ICI 2 

for 7^ |C|; < f ■ Finally, we calculate detHf. 



*IX/444(ICI/)/MICI/) 
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ICI? 



V ICI/ 



+ s(ICI/) 



444W* Vdci/) 
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lCl? S m 2 (|C|/ 



;i-|C|/cot(|C|z)) >0 



1/ 0111 y\s\D 

for 7^ |C|z < w. We conclude that the Hessian is positive definite and Re(7(z,w)) has a local 
minimum at w = iQ. Thus 

n 1 1 2 

Re(7(^^))>Re(7(x,^))U= iC = ^^-|C|/cot(|Cb)>c 2 |x| 2 if |C|/ < vr/4. 



l=i 



Put C = ez, then |C|z < ea and 1)7.11)1 holds with £0 = Jj- 

Estimate 1)7.12)1 is a consequence of estimates 1)7.10(1 and 1)7.11)1 since 



f(x, z,t + iez) = j(x, z,t + iez) -i ( T m + ie^)z m 

m=l 

3 

= j(x,z,t + iez) +e\z\ -i T m z m . 

m=l 
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□ 

Theorem 17.61 in a non-diagonal situation was proved in [I]. 
Lemma 7.7. If x is a non-zero vector in M. 4n , the integral 
(7-17) G(x,z)= I —^±—dr 



f 2n + 2 ( X ,Z,T) 

is absolutely convergent and one has for x ^ 

A G(x,z) = 0. 

Proof. Since the function V(r) does not depend on x and z, we have AqV = 0, X^V = 0, 
k = 1, 2, 3, 4, Z = 1, . . . , n, and the equation (|7.1|) reduced to the following one 

(7.18) Ao(V7" p ) = ~p(r p ~ VAq/ + (-p - l)Vf-v- 2 H(Xf)) . 

Here H(Xf) = \Xf\ 2 by Proposition 17.21 Moreover, taking into account that the complex 
action function f(x, z, r) satisfies the Hamilton-Jacobi equation (|6,3|) and p = 2n + 2, we get 
(7.19) 

A (Vr 2n - 2 ) = (-2n-2)(r 2n - 3 VA f + (-2n-3)Vf- 2n ~ 3 + (2n + 3)Vr 2n - 4 ^ r m ^- 

m=l 

Substituting the last term in the right hand side of (|7.19j) from the formula 

E /" (r m Vf- 2n - 3 ) = SVf- 2 - 3 + r 2n ~ 3 £ ^l^ 1 - (2n + 3)V J" 2 "" 4 £ r m |^, 

, OT m OT m ClT m 

m=l m=l m=l 

we deduce 

A Q (Vf- 2n - 2 ) = {-2n - 2)(/- 2 - 3 (y(A / - 2n) + £ r m |^) - £ o^-(r»V/-*- 8 )). 

m=l m=l 

Since the volume element V(r) is a solution of the transport equation (|7.7|) . finally, we obtain 
A G(x,z)= [ A (Vf- 2n - 2 )dT = (2n + 2) [ V (r m V/" 2n ~ 3 ) dr. 



We observe that 

(7.20) V(r) = TT — -► as one of the |r m | = i? -» oo, 

^ smh \t\i 

and 

n I ,2 I |2 

(7.21) |/|> ^M.| T | l00 th(|r|,)>H- 

because of \t\i coth(|r|/) > 1 for I = 1, . . . , n. The estimates (j7.2Uj) and (j7.21|) show 

(7.22) lim I J-(r mV f- 2n - 3 ) = 0. 

R -*°°J\T m \<R OT m K 

The last equality implies 

A (Vf- 2n ~ 2 ) = (2n + 2) / V ^- (r m y/- 2 ™- 3 ) dr = 0, 

^ 3 m=l aTm 
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that terminates the proof of Lemma 17.71 □ 

The argument of Lemma 17.71 is not valid for x = 0. In fact, the integral ()7.17|) is divergent 
because the denominator of the integrand contains — Yl m =i T mZ m which is zero along a hyper- 
plane of R 3 . We will treat the case by changing contour by adding a small imaginary part to 
the r m 's. We shall prove that for i/0we can change the contour in (|7.17|) and when x will 
be zero the integral (|7.17|) will still be convergent on the new contour. In order to achieve this 
goal, we need to use Proposition 17.61 

Proposition 7.8. For i / 0, the integral G(x, z) defined in \7.11\j is given by 
(7-23) G(x,z) = f ^j T + l£ ~ Z) . „ dr 



f 2n+2 (x, z,t + iez) 

for < e < £o sufficiently small. The integral makes sense even for x = and z / 0, 

so the function G(x,z) is well-defined (in fact, real analytic) except at the origin in R 4n x M 3 
and satisfies 

A G(x,z) = for (x,z)^(0,Q). 

Proof. We may prove this theorem by imitating the idea in [I]. Set 

Q, Ke = {£ = T + inz : rel 3 , \t\ <K, < n < e} C C 3 . 

Assume that |x| ^ 0. By Theorem 17.11 the differential form 

w = (V(()/f n+2 (x, z, C))d£i A #2 A dH 3 

is a homomorphic form of type (3, 0) in £Ir,e- It is easy to see that its differential is zero. 
Hence, by Stokes's Theorem 

/ oj = 0. 

J 8£Ik,e 

The boundary can be written as d£lx,e = <9^i U <9^2 U <9^3 • 

(1) The set dQi = {r € M 3 : \r\ < K, 7] = 0} which is such that 

lim [ ror,X^ \ dr = G ( x ' z ) 

since the integral (|7.17l) converges absolutely. 

(2) The set = {t + iez : |r| < K, r\ = e}. The integral (|7.23|) converges absolutely by 
Proposition 17.61 in this case. 

(3) The set d^lj, = {£ = r + iez : \t\ = K, < n < e}. Again, by Proposition 17.61 one has 

K^ooJlXh P n+2 (x,Z,T + l V z) 

By the discussion above, one may conclude that for i ^ 0, G(x, z) is given by the integral 
Q7.23JI on a shifted contour. Moreover, this integral is absolutely convergent even when x = 
and z 7^ by Proposition 17.61 We complete the proof of this theorem. □ 

Theorem 7.2. The kernel G(x,z) of the Green's function for the sub-Laplacian Aq is given 
by the formula 

, 2 2n {2iT) 2n + 3 f V{r + ie~z) 

G(x,zj = ; — - / -= , n . -dr. 

K ' (2n + l)! 7 R 3 / 2 «+ 2 (r + iez) 
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Proof. For any K > 0, denote 

G K (x,z) = — 1 / F(r + fez)dr f t^^e^dt. 

r(2n + 2) 7 R3 y 



The function Gk(x,z) is smooth everywhere and for 7^ (0,0), one has 

lir 

K-. 

Using Proposition 17.61 we know that 



lim Gk(x, z) = G(x, z). 

K—*oo 



\<cl 


air [ 


Jl 


13 JO 



sinh 2 |r|; 



2 



But the function (|x| 2 + |z|) 2n 2 £ L* oc (IR 4n+3 ) since the homogeneous degree is An + 6 in this 
case. It follows that 

lim G K (x,z) = G(x,z) in L} oc (R 4n+3 ) 

by the Dominated Convergence Theorem. We first calculate AoGr(x,z) for 1 / 0. We need 
to compute Ao(V^e - *^) which is 

AoiVe-V) = -tV • Ao(/)e"^ + t 2 H(x, z, Vf)Ve^. 

But the Hamilton-Jacobi equation ()6.3j) yields 

3 

A (Ve-tf) = - • A (/)e"^ + t 2 (/ - £ Tm ^-)Ve-^ 



m=l 



3 5 



= -te-tf(VAo(/) + 2 o— (^)) 

m=l 

m=l OTm 

We also know that 

[ R t 2n+3 e~^ fdt = -K 2n+3 e- K f + (2n + 3) /"* t 2n+2 e^dt. 
Jo Jo 



This implies that 

"A" 



(7.24) 



I dr [ t 2n+1 e^t 2 fVdt = -K 2n+3 [ Ve 
Jrz Jo Jr 3 



+ (2n + 3) t 2n+2 dt / e-^Vdr. 



Hence, 

-K P r 3 



A Q G K ( X ,z)=-—^— [ t 2 ^ 2 dt[ e-tf\j2^r n + (*of-2n-3)V 

^2n+3 

~ r(2n + 2 



l<-2n+3 

' Fe-^dr, 
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where we used (j7.22|) and ()7.24j) . The first integral vanishes since V satisfies the generalized 
transport equation (|7.7|) . Therefore, for i / 0, 

A G K (x,z) = - / V{r)e- K f^dr. 

T(2n + 2) J R3 

However, we may also change contour in the above integral and obtain 

jy-2n+3 f 

(7.25) A G K (x,z) = -——— V(T + iez)e- K ^ T+ie ^dT, for z + 0. 

T(2n + 2) 7 R3 

Since Gk(x,z) is smooth everywhere, the integral (|7.25j) provides the value of AqGk(x,z) at 
every point where the integral is convergent. It follows that AqGk(x, z) is equal to 



(2n + l)! 



everywhere and to — p|^pjyi Jr3 V(r)e K f( x > z > T )dT almost everywhere. Furthermore, for (x, z) 
in a compact set £/ with [/ disjoint from the origin, 

Re(/)(x, z, t + iez) > k > 0. 

Hence AqGk{x, z) — > uniformly as K ^ oo on compact subsets of M 4ri x M 3 which is 
disjoint from the origin. Now we need to compute the L 1 -norm of AqGk(x, z). Since the 
integral J K3 V {T)e~ t ^ x,z,T " > dr coincides with the integral of the right-hand side of (|7.25|) almost 
everywhere, we can just compute the following 

X:= _ R2n+3 ([[ VMe-KM^-^^^dTdzdx . 

(2n+l)! J R 4n J R 3 J R 3 

Since the integral converges absolutely, we may interchange the order of the integration by 
Fubini's Theorem. Let us integrate the z-variable first. 

-Kl 7 (x,T)-i-£ 3 m=1 TmZ m ] dz = f g-f Er=lNf[T[j«»th(|T|,)+«fE^=l'rm« n .] t i e 

Jr 3 



= - e -^V4Er=il-l?l-kcoth(H0 K -3 (27r) 3 Jc .-l (1)= ^-3 (27r) 3 (5(T) _ 

Here J r (<7)(^) = L3 g(x)e~ 2mx '^dx is the Fourier transform of the function g(x). We know that 

1 - 1 
lini7(x,r) = lim - ^ \ x \l \ T \i cot H\ T \l) = 4 1 



\x\ 2 
4 1 1 

1=1 



and 



It follows that 



rl 2 



lim V(t) = lim TT — ^ = 1. 

r^o r— >0 rr sinh |r|/ 



K- 3 (2vr) 3 / V(T)e- K ^ x '^5(T)d 
in the sense of distribution. Finally, one has 



A' 3 A' 2 " ' 
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This gives us 

K 2n+3 (2^)3 ^ (2vr) 2 "+3 

(2n + l)! K2n (2n + l)!' 

This proves that AoG\k(x, z) — > uniformly on compact sets on R 4n x R 3 disjoint from the 
origin with a constant integral over R 4n x R 3 . This means that when K — > oo, 



A Gjr(s,z) - -2 2 " o ' — J 



(2vr) 2 "+ 3 
(2n + l)!°(° m - 



On the other hand, Gk(x, z) — > G(x, 2) in L i 1 oc (R 4n+3 ) as if — > oo. Hence, 

A G^(x,z) -► A G(x,z) 
in the sense of distribution. Therefore, 

22n/-2 7r , )2n+3 

A G(,,,) = - 1 A_i_ r , (0)0) . 

The proof of the theorem is therefore complete. □ 

The symmetry of homogeneous H-type groups allows us to deduce another form of Green's 
function related to the homogeneous norm (see, for instance 

8. Estimates of the fundamental solution 

In this section, we discuss sharp estimates for the integral operator induced by the funda- 
mental solution G(x,z): 

G(g)(x,z) = G * g(x,z) = / G(y,w)g((y,w)~ 1 ■ (x, z))dydw 

in L p k (Q n ) Sobolev spaces, Hardy-Sobolev spaces H^(Q n ) for k £ Z + , < p < oo and Lipschitz 
spaces A^, Tp with (3 > 0. We consider here the sub-Laplacian Ao- It is easy to see from the 
group law that the operator is homogeneous of degree —2 under that non- isotropic dilation: 

5\ : (x,z) — ► (Xx,X 2 z). 

Hence, in general, the homogeneous degree of Q n is An + 6. As in Folland-Stein |14j and 
Koranyi j2U|, we may define a homogeneous norm 

n 4 3 i 

k^)i = ((EE4) 2 +£<r- 

3=1 k=l j=l 

Then we may define a pseudo metric by p(x, y) = |x • y" 1 ] where x = (x,z) and y = (y,w). 
Then it is easy to see that p is equivalent to the Carnot-Caratheodory metric. Using this 
metric, one may obtain estimates of G(x, z) in various function spaces. From the discussion 
in Section 7, we know that the fundamental solution G(x, z) with singularity at the origin is a 
homogeneous kernel of degree —An — 4. In fact, the fundamental solution G(x, z) satisfies the 
following size estimates: 



\G(x,z)\ < 



Ci|(x,z)| 2 8^G(x,z) ^ C 2 |(x,z)| 2 -H 



< 



and 



B((x,z),\(x,z)\)\' dx^---dxZ n ~ \B((x,z),\(x,z)\)y 
dh+b+toGfaz) < C 3 \(x,z)\ 2 ~ 2 ^+^+M 



dz% L dz% t d4 



B((x,z),\(x,z)\) 
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Therefore, it is a locally integrable function. Prom classical results, it is easy to see that the 
operator G originally defined on the Schwartz space S(Q n ) can be extended to a bounded 
operator from L p oc (Q n ) into L p oc (Q n ) for 1 < p < oo (see @] and ^1] ). Moreover, G is a 
smoothing operator. Hence, we may allow to differentiate the kernel G(x,z). Moreover, for 
any I = 1, . . . ,n and k = 1, ... ,4, X^iG originally defined on S(Q n ) can be extended to a 
bounded operator from L p oc (Q n ) into L p oc (Q n ) for 1 < p < oo. The problem now reduces to 
looking at the second derivatives (in the horizontal directions) of G(x, z). Before we go further, 
let us recall some basic definitions and properties of several functions spaces. 

• Lipschitz spaces As in we define the space Tp(Q n ) as the set of all bounded functions 
g with compact support on Q n such that 



sup 

(x,z),(j/,«;)eQ 71 



g((x,z)-(y,w)- 1 )-g(x,z) <C-\(y,w)f, < /3 < 1. 



When k < (3 < k + 1 with k = 1,2, we may define g E Tp(Q n ) as the set of all C k functions 
g with compact support on Q n such that V(X, X')g E Tg(Q n ), where V(X, X') is a monomial 
of degree k in vector fields X, X E V±. For integral value k, one may define the space Tg(Q n ) 
by interpolation. Furthermore, since \(x,z)\ < AWix^)]] 1 / 2 for \(x,z)\ small, one may conclude 
that Tg(Q n ) C A g / 2 (Q n )- Here ||(x,2i)|| is the Euclidean distance between the point (x,z) and 
the origin and A a (Q n ) is the isotropic Lipschitz space on Q n . As usual, the space A a (Q n ) is 
defined as the collection of all bounded functions g with compact support on Q n such that 



sup 

{x,z),(y,w)<=Q r ' 



9((x,z) ■ (y,w) - g(x,z) 



<(7.||(j/,«0H a , 0<a<l. 



For k < a < k + 1 with k = 1, 2, we may define g E A a (Q n ) is the set of all C k functions g 
with compact support on Q n such that V fc g E A a (Q n ). 

• Sobolev spaces. One may define the non-isotropic Sobolev spaces Su(Q n ) with k E Z + 
and 1 < p < oo as follows 

S P k (Q n ) = {/ : Q n - C : / E Lf(Q n ), V(X, X')f E L p (Q n )}, 

where V(X, X') is a monomial of degree k in vector fields X, X' E V\. Here L p (Q n ) is the L p 
Lebesgue space. 

• Hardy spaces. The Hardy space H p (Q n ) with < p < oo originally defined by maximal 
function as follows: a distribution / defined on Q n belongs H p (Q n ) if and only if the maximal 
function 

M(f)(x,z) = S up\f*^ £ \(x,z) E If(Q n ). 

e>0 

Here <f> E S(Q n ) with jg„ 0(x, z)dxdz = 1 and polyradial (see Chapter 4 in Folland-Stein fTS]). 
As usual, 

Mx,z) = e-( 4n+6 )0(e- 1 a;,e- 2 z). 
The space H p {Q n ) can be defined by atomic decomposition and maximal functions. 

Definition 8.1. A H p (Q n ) p-atom (0 < p < 1) is a compactly supported function a(x, z) such 
that the following conditions hold: 

(1) (size condition): there is a Q n -ball B Xo = B/ XOiZo \(r) = {x = (x,z) E Q n : p(xo,x) < r} 
whose closure contains supp(a) such that ||a(x)||x,oo < |i? XQ |~ 1 / p ; 

(2) (moment condition): 

a(x, z)V(x, z) dxdz = 



QUATERNION H-TYPE GROUPS Q n 37 



for all monomials V(x) = V(x, z) such that "P(x) = ( YYj=i x j 3 \ ' ' ' x jl ) z i z 2 z 3 anc ^ 



4 n 



r 1 

Y,^2 a ik + 2(/3i + 02 + ft) < [(4n + 6)(- - l) 
fc=i i=i p 

Here [s] is the integral part of s. 

Using the idea of atomic decomposition, we give the definition of H p (Q n ) as follows 

oo n 

H p (Q n ) = j/ e S'(Q n ) : / = ^ X k a k , where a fc are p-atoms, ^ |A fc | p < ooj, 

fc=l k=l 

and we define 

n 

n/iiW) = inf {Ei^i p }' 

/t=i 

where the infimum is taken over all possible atomic decompositions of /. Then "norm" 
||/||jyj>(Q«) is comparable to the t v norm of the sequence {X k } and the LP{Q n ) norm of the 
maximal function Ai(f). 

Now, following notations in [31, for g S S(Q n ), one has 

X jk X j > k >G{g)(x,z) = lim / (X jk G)(y,w)(X j , k ,g)((x,z)-(y,w)~ 1 )dydw 

e^0 + Jp(y,w)=e 

- lim / (X jk X jlk >G)(y,w)g((x,z) ■ (y^wy 1 ) dydw. 

£^0+ Jp(y,w)>£ 

The second term above is a generalized Calderon-Zygmund operator in the sense of |14j and 
Koranyi-Vagi |21j because Xj k Xj/ k >G is a kernel homogeneous of degree satisfying mean value 
zero property. An operator K is said to be a generalized Calderon-Zygmund operator if the 
following two conditions are satisfied: 

(i) K can be extended as a bounded operator on L 2 (Q n ), 

(ii) There is a sequence of positive constant numbers {Cj} such that for each j £ N, 
(8.1) 

, i 



_ 4»+li 



|K(x 1 ,x 2 )-K(x 1 ,x 3 )| 9 dx 1 <C r (2^(x 2 ,x 3 )) ~7 

2Jp(x2,x 3 )<p(xi,x 2 )<2J+ 1 p(x2,X3) / 



and 
(8.2) 



If \ 9 ■ 4n+6 

/ |K(x 2 ,x 1 )-ir(x3,x 1 )| 9 ( ix 1 <C i -(2V(x 2 ,x 3 )) "5^, 

\J2Jp(x2,X3)<p(xi,X2)<2J+ 1 p(x2,x 3 ) / 

here (q, q') is a fixed pair of positive numbers with | + ^- = 1 and 1 < </ < 2. We need the 
following theorem to complete our discussion on the estimates for G. The proof can be found 
in ^1] and |15j . 

Theorem 8.1. Let K be a generalized Calderon-Zygmund operator. Assume that the kernel 
satisfies conditions 118. 1\) - I8.fy) with {Cj} G I 1 . Then 

ll K (/)llLf(o n ) ^ c i»ll/llL»'(Q n )j 1 < p < °°; 

and 

{xGQ": |K(/)(x)| > A} 



C 



<-HfllLi(Q»), A>0. 
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When < p < 1, we may consider the boundedness of generalized Calderon-Zygmund 
operator acting on Hardy spaces H p (Q n ). We have the following theorem (see 

Theorem 8.2. Let K be a generalized Calderon-Zygmund operator. Suppose the following two 
conditions hold: 

(i) (kernel assumption): there exist s> [(4n + 6)(| — 1)] and e>^ — 1 such that 

{2 (4n+6)i(2 7 -l) (c ,. )2} G£ l wm 7=1 _I + £j 

p 

(ii) (adjoint operator assumption): K*(x"" • • • x^f z± 1 z§ 2 z% 3 ) = with 

££«M + £A<[(4n + 6)(i-l)], 

k=l 1=1 m=l P 

where K* is the adjoint operator of K . 

Then K can be extended as a bounded operator from H p (Q n ) into H p (Q n ) for < p < 1. 

We first make the following observation. For any polynomial V(X) of degree k in the 
horizontal vector fields X±i, . . . , X±i with I = 1, . . . , n, there exists another polynomial V(X) 
such that the following identity holds: 

(8.3) T(X)G = GP(X). 

Now we may apply the above theorems and (J8.3|) to conclude the following result. 

Theorem 8.3. The fundamental solution G(x, z) for the operator Ao defines an operator G 
which satisfies the following sharp estimates: 

i) XX'G defines a bounded operator from the Sobolev space S^.(Q n ) into Sobolev space 
S{{Q n ) for k E Z +; 1< p < oo and for all X, X' £ V i; 

ii) ZG defines a bounded operator from the non-isotropic Sobolev space S^{Q n ) into S^(Q n ) 
for k £ Z + , 1 < p < oo and for all Z £ V2; 

Hi) XX'G defines a bounded operator from H^,(Q n ) into H^(Q n ) for k G Z + and < p < 1 
and for all X, X' G V\; 

iv) ZG defines a bounded operator from H^(Q n ) into H^(Q n ) for k G Z + and < p < 1 
and for all Z G V2 ; 

v) G defines a bounded operator from Ap{Q n ) into Tp +2 (Q n ) n A / g + i(Q n ) for < (3 < 00. 
In conclusion, G gains two in horizontal directions and only gains one in missing directions. 
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